摘 要:为了解沙漠经济植物在新疆的潜在适生区,本文共选取了6种植物(板蓝根、甘草、黑果枸杞、红花、沙棘、枸杞)作为研究对象,基于6种植物在新疆区域的样本点以及生物气候、土壤和地形3类环境因子共31个环境变量,利用MaxEnt模型建立了6种沙漠经济植物的物种分布模型。在所建模型基础上,绘制了新疆6种沙漠经济植物在当前气候变化下的潜在适生区,并对6种沙漠经济植物在新疆的土地适宜性进行了可视化。最终构建了板蓝根、甘草、枸杞、黑果枸杞、红花、沙棘土地适宜性评价数据集6个。结果显示各植物适宜性评价模型的AUC值均大于0.8,表明预测结果较为准确。每个数据集均包含2个数据文件:(1) 该类植物土地适宜性评价结果数据,以.txt和.tif格式储存;(2) 该类植物空间分布点数据(引用数据),包括.xlsx和.shp格式。
关键词:MaxEnt模型;沙漠经济植物;空间分布;适生区
DOI: https://doi.org/10.3974/geodp.2025.02.09
CSTR: https://cstr.escience.org.cn/CSTR:20146.14.2025.02.09
数据可用性声明:本文关联实体数据集已在《全球变化数据仓储电子杂志(中英文)》出版,可获取:
https://doi.org/10.3974/geodb.2025.02.02.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2025.02.02.V1.
https://doi.org/10.3974/geodb.2025.02.03.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2025.02.03.V1.
https://doi.org/10.3974/geodb.2025.02.04.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2025.02.04.V1.
https://doi.org/10.3974/geodb.2025.02.05.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2025.02.05.V1.
https://doi.org/10.3974/geodb.2025.02.06.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2025.02.06.V1.
https://doi.org/10.3974/geodb.2025.02.07.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2025.02.07.V1.
物种适生区预测是生态学领域中的关键研究课题,它在多个方面发挥着重要作用,包括有害生物的防治[1]、珍稀濒危物种的保护[2,3]以及作物引种和栽培[4,5]等。物种分布模型(species distribution models,SDMs)是预测物种潜在地理分布和栖息地适生性的有利工具[6–8]。最大熵模型(MaxEnt)是物种分布模型之一,该模型是一种通用机器学习方法,利用已知的物种分布记录及其相关的环境变量来推断物种的生态需求,进而预测研究区域内物种的潜在适生分布区[9,10]。因MaxEnt模型具有样本需求量低、操作简便、预测精度高等特点,而被广泛认可和使用[11]。例如:张亮[12]利用MaxEnt模型对黑果枸杞(Lycium ruthenicum Murray)在中国的潜在适生区进行了预测;陈松清[13]基于MaxEnt模型对新疆沙棘(Hippohgae rhamnoides L.)的地理分布做了预测。此外,该模型还被用于预测农作物的潜在分布区域[14]、评估气候变化对物种生长的影响[15]、以及优化物种种植结构布局[16]等。
目前,板蓝根(Isatis tinctoria L.)、红花(Carthamus tinctorius L.)、枸杞(Lycium dasystemum Pojark.)的适生区研究较少,且缺乏新疆区域适生区的可视化分析。甘草(Glycyrrhiza uralensis Fisch.)、黑果枸杞、沙棘的适生区预测主要集中在全中国范围内,新疆区域内适生区的划分仍可一步细化[12,17,18]。综合来看,该6种沙漠经济植物在新疆的适生区预测仍需进一步研究。
本文以板蓝根、甘草、枸杞、黑果枸杞、红花和沙棘为研究对象,基于MaxEnt模型和ArcGIS技术,利用全球生物多样性信息网络数据库(Global Biodiversity Information Facility,GBIF)、WorldClim2.1数据库和世界土壤数据库获取植物样本点数据和环境变量数据,分析了6种沙漠经济植物适生区在新疆的空间分布,以期为这些植物在新疆的合理利用和规划提供参考。
新疆板蓝根土地适宜性评价数据集[19]、新疆甘草土地适宜性评价数据集[20]、新疆枸杞土地适宜性评价数据集[21]、新疆黑果枸杞土地适宜性评价数据集[22]、新疆红花土地适宜性评价数据集[23]、新疆沙棘土地适宜性评价数据集[24]的名称、作者信息、地理区域、数据年代、数据集组成、数据出版与共享服务平台、数据共享政策等信息见表1。
| 条目 | 描述 |
|---|---|
| 数据集名称 | 新疆板蓝根土地适宜性评价数据集 |
| 数据集短名 | Isatis tinctoria L._XJSuitable |
| 作者信息 | 张萍,中国科学院新疆生态与地理研究所,z1571824849@163.com 俞祥祥,中国科学院新疆生态与地理研究所,yuxiangxiang@ms.xjb.ac.cn 常存,中国科学院新疆生态与地理研究所,changcun@ms.xjb.ac.cn 张恒,中国科学院新疆生态与地理研究所,zhangheng@ms.xjb.ac.cn 范敬龙,中国科学院新疆生态与地理研究所,fanjl@ms.xjb.ac.cn |
| 地理区域 | 新疆 |
| 数据年代 | 生物气候变量1970–2000,地形数据2000,土壤数据2009 |
| 时间分辨率 | 年 |
| 空间分辨率 | 30″ |
| 数据格式 | .shp、.tif、.xlsx、.txt |
| 数据量 | 292 KB |
| 数据集组成 | 板蓝根土地适宜性评价结果数据,空间分布点数据(引用数据) |
| 数据集名称 | 新疆甘草土地适宜性评价数据集 |
| 数据集短名 | Glycyrrhiza uralensis Fisch._XJSuitable |
| 作者信息 | 张萍,中国科学院新疆生态与地理研究所,z1571824849@163.com 俞祥祥,中国科学院新疆生态与地理研究所,yuxiangxiang@ms.xjb.ac.cn 常存,中国科学院新疆生态与地理研究所,changcun@ms.xjb.ac.cn 马强,中国石油天然气股份有限公司塔里木油田分公司,462920056@qq.com 王龙,中国石油天然气股份有限公司塔里木油田分公司,41063307@qq.com 范敬龙,中国科学院新疆生态与地理研究所,fanjl@ms.xjb.ac.cn |
| 地理区域 | 新疆 |
| 数据年代 | 生物气候变量1970–2000,地形数据2000,土壤数据2009 |
| 时间分辨率 | 年 |
| 空间分辨率 | 30″ |
| 数据格式 | .shp、.tif、.xlsx、.txt |
| 数据量 | 313 KB |
| 数据集组成 | 甘草土地适宜性评价结果数据,空间分布点数据(引用数据) |
| 数据集名称 | 新疆枸杞土地适宜性评价数据集 |
| 数据集短名 | Lycium dasystemum Pojark._XJSuitable |
| 作者信息 | 张萍,中国科学院新疆生态与地理研究所,z1571824849@163.com 俞祥祥,中国科学院新疆生态与地理研究所,yuxiangxiang@ms.xjb.ac.cn 常存,中国科学院新疆生态与地理研究所,changcun@ms.xjb.ac.cn 张恒,中国科学院新疆生态与地理研究所,zhangheng@ms.xjb.ac.cn 范敬龙,中国科学院新疆生态与地理研究所,fanjl@ms.xjb.ac.cn |
| 地理区域 | 新疆 |
| 数据年代 | 生物气候变量1970–2000,地形数据2000,土壤数据2009 |
| 时间分辨率 | 年 |
| 空间分辨率 | 30″ |
| 数据格式 | .shp、.tif、.xlsx、.txt |
| 数据量 | 285 KB |
| 数据集组成 | 枸杞土地适宜性评价结果数据,空间分布点数据(引用数据) |
| 数据集名称 | 新疆黑果枸杞土地适宜性评价数据集 |
| 数据集短名 | Lycium ruthenicum Murray_XJSuitable |
| 作者信息 | 张萍,中国科学院新疆生态与地理研究所,z1571824849@163.com 俞祥祥,中国科学院新疆生态与地理研究所,yuxiangxiang@ms.xjb.ac.cn 常存,中国科学院新疆生态与地理研究所,changcun@ms.xjb.ac.cn 张恒,中国科学院新疆生态与地理研究所,zhangheng@ms.xjb.ac.cn 范敬龙,中国科学院新疆生态与地理研究所,fanjl@ms.xjb.ac.cn |
| 地理区域 | 新疆 |
| 数据年代 | 生物气候变量1970–2000,地形数据2000,土壤数据2009 |
| 时间分辨率 | 年 |
| 空间分辨率 | 30″ |
| 数据格式 | .shp、.tif、.xlsx、.txt |
| 数据量 | 332 KB |
| 数据集组成 | 黑果枸杞土地适宜性评价结果数据,空间分布点数据(引用数据) |
| 数据集名称 | 新疆红花土地适宜性评价数据集 |
| 数据集短名 | Carthamus tinctorius L._XJSuitable |
| 作者信息 | 张萍,中国科学院新疆生态与地理研究所,z1571824849@163.com 俞祥祥,中国科学院新疆生态与地理研究所,yuxiangxiang@ms.xjb.ac.cn 常存,中国科学院新疆生态与地理研究所,changcun@ms.xjb.ac.cn 张恒,中国科学院新疆生态与地理研究所,zhangheng@ms.xjb.ac.cn 范敬龙,中国科学院新疆生态与地理研究所,fanjl@ms.xjb.ac.cn |
| 地理区域 | 新疆 |
| 数据年代 | 生物气候变量1970–2000,地形数据2000,土壤数据2009 |
| 时间分辨率 | 年 |
| 空间分辨率 | 30″ |
| 数据格式 | .shp、.tif、.xlsx、.txt |
| 数据量 | 231 KB |
| 数据集组成 | 红花土地适宜性评价结果数据,空间分布点数据(引用数据) |
| 数据集名称 | 新疆沙棘土地适宜性评价数据集 |
| 数据集短名 | Hippophae rhamnoides L._XJSuitable |
| 作者信息 | 张萍,中国科学院新疆生态与地理研究所,z1571824849@163.com 俞祥祥,中国科学院新疆生态与地理研究所,yuxiangxiang@ms.xjb.ac.cn 常存,中国科学院新疆生态与地理研究所,changcun@ms.xjb.ac.cn 张恒,中国科学院新疆生态与地理研究所,zhangheng@ms.xjb.ac.cn 范敬龙,中国科学院新疆生态与地理研究所,fanjl@ms.xjb.ac.cn |
| 地理区域 | 新疆 |
| 数据年代 | 生物气候变量1970–2000,地形数据2000,土壤数据2009 |
| 时间分辨率 | 年 |
| 空间分辨率 | 30″ |
| 数据格式 | .shp、.tif、.xlsx、.txt |
| 数据量 | 378 KB |
| 数据集组成 | 沙棘土地适宜性评价结果数据,空间分布点数据(引用数据) |
| 基金项目 | 新疆维吾尔自治区(2022B03030) |
| 数据计算环境 | MaxEnt,ArcGIS |
| 出版与共享服务平台 | 全球变化科学研究数据出版系统 http://www.geodoi.ac.cn |
| 地址 | 北京市朝阳区大屯路甲11号100101,中国科学院地理科学与资源研究所 |
| 数据共享政策 | (1)"数据"以最便利的方式通过互联网系统免费向全社会开放,用户免费浏览、免费下载;(2)最终用户使用"数据"需要按照引用格式在参考文献或适当的位置标注数据来源;(3)增值服务用户或以任何形式散发和传播(包括通过计算机服务器)"数据"的用户需要与《全球变化数据学报(中英文)》编辑部签署书面协议,获得许可;(4)摘取"数据"中的部分记录创作新数据的作者需要遵循10%引用原则,即从本数据集中摘取的数据记录少于新数据集总记录量的10%,同时需要对摘取的数据记录标注数据来源[25]。 |
| 数据和论文检索系统 | DOI,CSTR,Crossref,DCI,CSCD,CNKI,SciEngine,WDS,GEOSS,PubScholar,CKRSC |
研究所选用的植物分布数据来源于全球生物多样性信息网络数据库(GBIF)[1],采样年限为1949–2020。基于此平台,本研究获取新疆范围内的板蓝根、甘草、黑果枸杞、红花、沙棘和枸杞采样分布点[26]的精确经纬度信息,作为各类植物的空间分布基础数据。
研究共选取31个环境数据用于建模,包括19个生物气候变量(Bio1–Bio19)、9个土壤变量以及海拔、坡度和坡向变量。其中,19个生物气候变量及海拔数据全部来自WorldClim2.1数据库[27];坡度坡向数据基于海拔数据利用ArcGIS计算得到;表层土壤数据来源于国家冰川冻土沙漠科学数据中心[28]。具体环境变量见表2。
| 类型 | 变量 | 描述 | 单位 |
|---|---|---|---|
| 生物气候变量 | Bio1 | 年平均气温 | ℃ |
| Bio2 | 昼夜温差月均值 | ℃ | |
| Bio3 | 等温性((Bio2/Bio7)×100) | \ | |
| Bio4 | 气温季节变化标准差(标准差×100) | \ | |
| Bio5 | 最暖月份最高温度 | ℃ | |
| Bio6 | 最冷月份最低温度 | ℃ | |
| Bio7 | 气温年较差(bio5–bio6) | ℃ | |
| Bio8 | 最湿季度平均温度 | ℃ | |
| Bio9 | 最干季度平均温度 | ℃ | |
| Bio10 | 最暖季度平均温度 | ℃ | |
| Bio11 | 最冷季度平均温度 | ℃ | |
| Bio12 | 年降水量 | mm | |
| Bio13 | 最湿月份降水量 | mm | |
| Bio14 | 最干月份降水量 | mm | |
| Bio15 | 降水量季节性变化 | \ | |
| Bio16 | 最湿季度降水量 | mm | |
| Bio17 | 最干季度降水量 | mm | |
| Bio18 | 最暖季度降水量 | mm | |
| Bio19 | 最冷季度降水量 | mm | |
| 地形 | Elev | 海拔 | m |
| Slope | 坡度 | ° | |
| Aspect | 坡向 | \ | |
| 土壤特性 | T_CaCO3 | 表层土壤碳酸钙(CaCO3)含量 | % |
| T_cec_soil | 表层土壤阳离子交换容量 | cmol/kg | |
| T_clay | 表层土壤粘粒含量 | % | |
| T_esp | 表层土壤酸碱度 | \ | |
| T_gravel | 表层土壤砾石含量 | % | |
| TOC | 表层土壤有机碳含量 | % | |
| T_pH_H2O | 表层土壤pH值 | –log(H+) | |
| T_sand | 表层土壤砂粒含量 | % | |
| T_texture | 表层土壤质地 |
本研究中,沙漠经济植物的空间范围为新疆。环境因子图层的空间分辨率为30弧秒,同一栅格内各处的气候因子数据相同,同一栅格内的冗余分布数据或重复分布数据会导致模型的过度拟合。因此,选用ENMTools工具对所收集的分布数据进行筛选。ENMTools工具是一款用于物种分布建模和预测的开源软件,它提供了一系列功能强大的工具,可以帮助研究人员从大量的物种分布数据中提取有用的信息。该工具可以自动匹配分析所用的环境因子图层分辨率大小,删除同一栅格内的冗余数据,使模型拟合结果更好。
环境因子之间存在相关性,全部用于模型预测会导致过度拟合。正确量化已有物种分布点与对应环境变量间的关系是物种分布模拟的基础。模拟物种分布时应选择对物种分布有重要影响的环境变量同时避免环境变量较高共线性导致模型预测结果过拟合[29]。因此,为避免环境变量间共线性导致过拟合,使用ENMTools软件进行环境变量的相关性分析,得到各因子之间的相关性系数(图1)。相关性系数是一种定量衡量变量之间的相关系数,其大小在−1和1之间,反映了两个变量之间线性相关的程度,‒1为完美的负相关,1为完美的正相关,绝对值越趋于1,则说明变量之间的相关性越强。将全部31个环境变量进行建模,得到各环境变量的贡献率,通过每个环境变量的贡献率和相关性系数[30],进一步剔除相关系数的绝对值≥0.8且贡献率<0.5的环境因子,使用剩余的环境变量进行最终建模。
图1 环境因子相关性分析结果图
下载将处理后的物种分布数据和环境因子数据导入MaxEnt模型中进行训练和测试。对于每个物种,使用交叉验证提取测试集。将所有的发生记录分为10个子集,其中1个子集作为测试集,其余9个子集作为训练集,设置最大迭代次数为10,000,数据输出格式为logistic。选择刀切法预测变量的重要性,并生成响应曲线。通过10次重复的交叉验证来建立模型,以减少不确定性。取10次结果的平均值作为最终结果。
模型的精度检验主要采取MaxEnt模型输出结果的受试者工作特征曲线(ROC)下的面积AUC值进行评估,ROC曲线显示了在不同阈值下正确预测的分布点与错误的预测点之间的关系。AUC的取值范围为0‒1,数值越大表示预测结果的可信度越高,AUC值低于0.8表示可信度低,0.8‒0.9表示准确度好,0.9‒1.0表示准确度极好[31]。
在ArcGIS 10.8软件中加载MaxEnt模型运算结果,利用空间分析模块中的重分类功能对生成的栅格数据文件进行重分类,绘制6种沙漠经济植物在新疆范围的潜在适生区分布图。MaxEnt模型预测的适生值是连续的栅格数据,数值在0‒1之间。植物适宜性等级划分标准按自然断点法[32]划分为不适生区、低适生区、中适生区、高适生区。每种植物适生值自然断点法划分标准如表3所示。
| 植物名称 | 不适生区 | 低适生区 | 中适生区 | 高适生区 |
|---|---|---|---|---|
| 板蓝根 | 0‒0.06 | 0.06‒0.23 | 0.23‒0.49 | 0.49‒1 |
| 甘草 | 0‒0.09 | 0.09‒0.26 | 0.26‒0.51 | 0.51‒1 |
| 黑果枸杞 | 0‒0.08 | 0.08‒0.25 | 0.25‒0.52 | 0.52‒1 |
| 红花 | 0‒0.07 | 0.07‒0.22 | 0.22‒0.49 | 0.49‒1 |
| 沙棘 | 0‒0.09 | 0.09‒0.26 | 0.26‒0.54 | 0.54‒1 |
| 枸杞 | 0‒0.07 | 0.07‒0.23 | 0.23‒0.50 | 0.50‒1 |
综上,本研究技术路线图如图2所示。
图2 数据集研发技术路线图
下载共计6个数据集,每个数据集包含2个数据文件:(1) 新疆沙漠经济植物(板蓝根、甘草、黑果枸杞、红花、沙棘、枸杞)土地适宜性评价结果数据,以.txt和.tif格式储存;(2) 新疆沙漠经济植物(板蓝根、甘草、黑果枸杞、红花、沙棘、枸杞)空间分布点数据(引用数据),包括.xlsx和.shp格式。
经用ENMTools工具对所收集的分布数据进行筛选,最终选取板蓝根、甘草、黑果枸杞、红花、沙棘和枸杞植物样本点数量分别为22、36、43、19、66和16。根据环境变量相关性分析结果和贡献度对环境变量进行过滤,过滤后板蓝根、甘草、黑果枸杞、红花、沙棘和枸杞的环境变量分别剩余13、15、16、13、15和15个,如图3各分图所示。板蓝根、甘草、黑果枸杞、红花、沙棘和枸杞适宜性评价模型的AUC值分别为:0.950,0.874,0.877,0.905,0.896,0.943。每种植物的AUC值均大于0.8,MaxEnt模型预测结果良好,且板蓝根、红花和枸杞的AUC值均大于0.9,表明预测结果较准确。
图3 6种沙漠经济植物模型预测结果AUC值对比图
下载当前气候条件下,板蓝根、甘草、黑果枸杞、红花、沙棘和枸杞在新疆的潜在适生区如图4所示,6种沙漠经济植物在新疆的适生区主要分布在绿洲区域,但各植物的潜在分布区存在一定的差异。具体如下:
图4 新疆6种沙漠经济植物的土地适宜性评价图
(该图基于新疆维吾尔自治区自然资源厅标准地图服务网站新S(2023)061号的标准地图制作,底图边界无修改)
(1)板蓝根高适生区主要位于新疆的西北部,集中在伊犁州、塔城地区和阿勒泰地区;中适生区分布在塔城地区、伊犁州和昌吉州等地带;低适生区主要位于哈密地区、南疆部分绿洲区和古尔班通古特沙漠等地;不适生区大多分布在南疆(图4a)。
(2)甘草高适生区集中在北疆,呈片状分布;中适生区分布在高适生区外围;低适生区主要分布在哈密地区、巴州、阿克苏和喀什等地;不适生区集中在南疆,位于巴州南部以及和田地区(图4b)。
(3)黑果枸杞高适生区呈现南疆高于北疆的分布特点,且呈带状分布,主要分布在南疆的和田-喀什-克州-阿克苏一带的绿洲区,北疆则主要分布在博州-塔城-昌吉州-乌鲁木齐-吐鲁番一带的绿洲带;中适生区主要集中在塔城地区;低适生区主要分布在巴州和阿勒泰东部(图4c)。
(4)红花高适生区主要位于各地州绿洲区,但分布稀疏且零散;中适生区集中在塔城地区的西北部和阿勒泰地区西部;低适生区主要分布在塔城地区;不适生区的范围较广,主要分布在新疆东部,特别是在阿勒泰、哈密、巴州和和田等地(图4d)。
(5)沙棘高适生区在各地州绿洲区都有分布,且呈片状分布,主要位于伊犁州和喀什地区;中适生区主要分布在塔城地区;低适生区主要分布在哈密地区、克州和喀什南部;不适生区大多集中在南疆,主要在巴州南部及和田北部(图4e)。
(6)枸杞高适生区主要位于北疆,基本呈片状分布,南疆的和田和喀什地区亦有一定分布;中适生区主要分布在塔城地区和阿勒泰地区;低适生区主要集中在和田西部、喀什东部、阿勒泰地区西部和哈密地区;不适生区主要在巴州南部(图4f)。
综合来看,在新疆地区,6种沙漠经济植物的高适生区和中适生区范围普遍较小,其中高适生区的占比均不超过5%。高适生区占比较高的植物为甘草、沙棘和枸杞,分别为3.61%、2.58%和2.46%。大部分区域为不适生区,其中红花的不适生区占比最高,达到86.27%;其次是板蓝根和黑果枸杞,不适生区占比分别为84.80%和75.01%;甘草的不适生区占比最低,为69.91%。从各类植物的低适生区范围来看,占比最高的植物是枸杞,为17.98%,最低的植物是红花,占比为9.39%。
本研究基于MaxEnt模型和ArcGIS,结合生物气候变量、地形和土壤多维环境变量,构建了新疆沙漠经济植物的适生区预测模型。通过自然断点法将适生区划分为4个等级,所有植物适宜性评价模型的AUC值均超过0.8,验证了模型预测结果的生态合理性。研究结果表明,新疆6种沙漠经济植物的潜在适生区呈现明显的空间异质性,其分布格局与绿洲生态系统密切关联。高适生区占比普遍低于5%,其中甘草(3.61%)、沙棘(2.58%)和枸杞(2.46%)具有相对优势。
本研究构建的新疆沙漠经济植物适生区数据集,可为沙漠产业空间规划提供参考,为植物资源的可持续利用提供了科学依据。未来的研究可进一步探讨沙漠经济植物适生区在未来气候变化下的时空演变过程,并量化未来各植物适生区质心迁移速率,以评估未来沙漠植物资源的可持续利用潜力。
作者分工:张萍收集和处理了数据并撰写了论文,王欣鹏、张冠峰和张恒参与处理了数据,范敬龙对数据集开发做了总体设计,张恒和俞祥祥指导了本文的写作。
利益冲突声明:本研究不存在研究者以及与公开研究成果有关的利益冲突。
[1] GBIF. http://www,gbif.org.